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Abstract 

It is possible to formulate fluid dynamics in terms of group-valued variables. This is particu- 
larly suited to the cases where the fluid has nonabelian charges and is coupled to nonabelian 
gauge fields. We explore this formulation further in this paper. An action for a fluid of 
relativistic particles (with and without spin) is given in terms of the Lorentz and Poincare 
(or de Sitter) groups. Considering the case of particles with flavor symmetries, a general 
fluid action which also incorporates all flavor anomalies is given. The chiral magnetic and 
chiral vorticity effects as well as the consequences of the mixed gauge-gravity anomaly are 
discussed. 



1 Introduction 



The description of fluid dynamics, especially for systems made of particles carrying non- 
abelian charges, has become an important research topic with the discovery of of the state 
of unconfined quark and gluons, the quark-gluon plasma, in heavy ion collisions. Field the- 
oretic analyses, augmented with Boltzmann-type kinetic equations, can be used to "derive" 
the equations of fluid dynamics, but are generally limited to dilute systems near equilib- 
rium. However, the basic equations can be formulated using general principles and therefore 
have a regime of validity significantly beyond the context of the derivation based on kinetic 
equations. The question of a derivation based on symmetry principles generalizing the usual 
equations of magnetohydrodynamics to include nonabelian charges and fields is interesting 
in its own right, but has enhanced relevance after the discovery of the quark-gluon plasma. 
Such an approach was developed in [U [2], where the fluid degrees of freedom were shown 
to be naturally described by the elements of a Lie group. The method applies to ordi- 
nary hydrodynamics as well, but becomes particularly useful for incorporating nonabelian 
symmetries and coupling to nonabelian fields. 

A new impetus to such analyses has come from the recent work on the chiral magnetic 
effect [3]. The specific case of interest has been the charge separation and a corresponding 
electromagnetic current induced by the axial anomaly, which can be demonstrated by the 
standard diagrammatic techniques. An interesting question to ask is then: Is there an 
effective description of the anomalies and how they affect the fluid dynamics? A related 
question is one of generalization to all flavor anomalies, even though they may not be 
of immediate relevance to the quark-gluon plasma. Symmetries are obviously front-and- 
center in analyzing anomalies and so our approach to fluid dynamics based on group-valued 
variables would seem tailor-made for these questions. This is exactly the subject of the 
present paper. 

The chiral magnetic effect, we may note, has led to a significant body of literature on 
related topics. The possibility of describing the chiral magnetic effect using hydrodynamics 
and thermodynamics is explored in Transport in arbitrary dimensions induced by 
anomalies has also been discussed in 0. Since there is considerable evidence that the quark- 
gluon plasma can be described as a strongly coupled fluid, the holographic correspondence 
can provide another method towards its analysis. The description of the chiral magnetic 
effect using holographic approaches such as the AdS/CFT correspondence or the Sakai- 
Sugimoto model is given in There is also an attempt to understand the chiral magnetic 
effect [7\ using the fluid/gravity correspondence of [8j. For completeness, we also note that 
there have been many lattice simulations of the chiral magnetic effect [S]. 

The focus in most of the literature has been on the computation of transport coefficients 
or the modifications of the energy- momentum tensor and the currents, and then the sub- 
sequent incorporation of these in the equations of motion of fluid dynamics. Our approach 
will be to write an effective action for anomalies directly in the fluid language, in other 
words, we obtain the fluid version of the Wess-Zumino term for anomalies. The action- 
based approach gives a simple starting point for all flavor anomalies. We also discuss some 
aspects of the mixed gauge-gravity anomaly in the standard model. (The mixed anomalies. 
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and the possibility of a chiral vortex effect, are also considered from the point of view of 
Kubo formulas and transport coefficients in |10).) On the negative side, the action-based 
approach will not include dissipative effects; they have to be added on after the equations 
of motion have been obtained by the variational principle. We note that an effective ac- 
tion approach has been given in two dimensions [llj, although the formulation is still very 
different from ours. 

In section 2, we give a brief resume of the formulation of ordinary, nonrelativistic or 
relativistic, fluid dynamics in terms of group variables. We then describe how nonabelian 
internal symmetries are included to obtain a nonabelian magnetohydrodynamics. In section 
3, we follow a similar approach to construct the action for a fluid of spinning particles in 
terms of the Lorentz group or Poincare (or de Sitter) groups, the latter being adaptable 
to the spinless case as well. The fluid description for the quarks in the standard model is 
given in section 4, taking a fluid of the up, down and strange quarks as an example. The 
full fluid action for these degrees of freedom including anomalies is given in this section. 
The standard chiral magnetic effect, the chiral vorticity effect and mixed gauge- gravity 
anomalies are discussed in this framework in section 5. 



2 Lagrangians and perfect fluids: a short resume 



We start with a recapitulation of the formulation of hydrodynamics in terms of group theory. 
We will be brief, since this is reviewed in detail in [T| l2]. Ordinary fluid dynamics can be 
viewed as a Poisson bracket system with 



6p * \6vi 



—d- (— 

6p * \6vi 



for F, G which are functions of the density p and fluid velocity Vi. The Hamiltonian, 

'1 



H 



■p^' + V{p) 



(1) 



(2) 



is then easily verified, via the brackets ([T]) to lead to the continuity and Euler equations, 
the fluid pressure being P = p^ — V. 

The difficulty with this framework is also well known. The helicity C, defined by, 

C=— I e'^'^vidjvk (3) 



1 

8^ 



is seen to Poisson commute with all observables, i.e. [F, C] = for all F. Viewing this 
from a quantum point of view, we see that the values of C are superselected. It is therefore 
necessary to specify a value for C and consider the restricted Hamiltonian dynamics for that 
sector itself. Alternatively, if we think of the Poisson brackets to be written as [^aj Cfe] = ^aby 
for ^ being p and Vi„ then the symplectic structure is obtained as the inverse of Kab- (Usually 
in starting from a Lagrangian, we obtain the symplectic structure and invert it to obtain the 
Poisson brackets.). Since C commutes with all observables, we see that {6C/6vi) is a zero 
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mode for K and hence we cannot relate Eqs. (|ll2p to a symplectic structure or Lagrangian 
description without first restricting the value of C. Thus, to obtain a Lagrangian description 
we must first fix C and then seek a parametrization for Vi which does not further change 
the value of C. This is given by the Clebsch parametrization, 



Vi = di9 + a difi 



(4) 



for an arbitrary functions 9,a, /3. In this case C = if 0, a, /3 are single- valued nonsingular 
functions vanishing at infinity. (They are known as the Monge potentials.) A suitable action 
for fluid dynamics is then |2l [T2] 



p9 + pa j3 



1 2 
- p V 

2 



V 



(5) 



We note that {p,0), {pa, (5) form two sets of canonically conjugate pairs. 



Now we introduce an element g of the group SU{1, 1) which may be parametrized in 
general as. 



1 



1 u 



e-'^^l'^ 



- uu 1 

where « is a complex variable. Direct calculation shows that 



a 



Tr (ct3 dg) = dO + a d/3, 
/3=H/2)ln(|). 



(6) 



(7) 



uu 



1 



uu 



The ^-direction in g corresponds to the compact direction, the U{1) subgroup generated by 
|(73, while a and /3 parametrize SU{1, 1)/U{1). The action can now be written as. 



S 



Jiji 



+ V 



(8) 



where we denote = p. The elimination of jj in ([8|) leads to the version ([5]). 

The relativistic generalization of fluid dynamics and the action ([8]) is also very straight- 
forward. It is given by 

S = -i j fl^ {asg-^d^g) - F{n) (9) 
where Fin) is a function of the variable n, which is defined by j^j^j, = n?- Equivalently 



(10) 



where n'^ is a four- vector obeying n'^n^ = 1. It may be considered as the four- velocity of 
the fluid and n identifled as the invariant density. The energy-momentum tensor for ([9]) has 
the perfect fluid form 

T^"' = nF'u^'u" - g^'^nF' - F) (11) 
identifying the pressure as P = nF' — F. The function F is thus the enthalpy. 



4 



We have obtained a group theoretic formulation of ordinary fluid dynamics. The exis- 
tence of a compact direction, namely, the U{1) direction of the SU{1, 1), may seem a little 
puzzling at first, since at the level of the classical Clebsch parametrization, this was not a 
requirement. The Poisson bracket obtained from ([8]) or ([9|) gives 

[pif),g{x)] = -ig{x)^f{x). (12) 

This means that in the quantum theory 

gU = ge'^''^ = -g (13) 

for U = exp [— 27ri J p] . Since all observables involve even powers of g, they are invariant 
under the action of U. This means that we can set U = 1, giving J p = N for some integer 
A^. The existence of the compact direction thus requires the quantization of / p in the 
quantum theory; this is equivalent to saying that the fluid is made of particles with p being 
the particle density [13J. Thus, rather than a defect of the group-theoretic parametrization 
([7]) in comparison to the classical Clebsch parametrization (jj]), we view this as a good 
feature of the description in ([8]), Q. [If vorticity were also quantized we would use SU{2) 
in place of SU (1, 1).] 

It is now easy enough to obtain the generalization to carrying nonabelian gauge charges, 
corresponding to a compact Lie group G. 

First consider SU{2). At the particle level, the dynamics of a particle carrying SU{2) 
charges is given by the Wong equations which have the action j l4t IT5] 



S 



- m + A^Q^Xi - i w Tr(cJ3C/~^ g) 



(14) 



where = Tr{g(T^g-^e), = ^a". 



The last term in (|14p is the co-adjoint orbit action which describes the dynamics of 
the gauge charges and which, upon quantization, gives the Hilbert space corresponding to 
one unitary irreducible representation (UIR) of SU{2) corresponding to the highest weight 
w/2, hence of dimension w + 1. Q"^ then become operators realizing the charge algebra 

Q",(3^] =ir^'Q'- (15) 

Under g — )• gexp {ia3(j)/2), the change in the action is given by AS = wAcj). Thus single- 
valuedness of e*'^ when cj) traces out a closed path in SU{2) leads to the quantization of w. 
The crucial co-adjoint orbit term, when generalized to several particles, becomes 



S 



I J dtY,wxTi{a3g^'gx) (16) 



where we have a separate g for each A, and likewise for w, with A indexing the particles. 
The continuum limit of (jl6p may be taken, as one does for the Lagrange approach to fluids, 
by A — )■ X, — J d^:x./v, w\/v — >• /9(x). This leads to 



S 



\ j <fx pTria^g-^g) (17) 
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where g = g{x,t). Taking this as the leading term, namely, as the term responsible for the 
symplectic structure, we can write an action 

S = -iJ d^x f Tr {asg'^D^^g) - J F{n) + Sym (18) 

where D^g = d^g + A^g, = -if^A^^, T = aV2. 

The velocity for the transport of the nonabelian charge can be introduced via = nu^, 
ti^ = 1 . The current which couples to the SU (2) gauge field A'^ is given by 

= Tr (^3 g-^f'g) f = Q'^f (19) 

which is in the Eckart form [16] . Starting with the action, one can easily verify the following 

mm- 

1. The equations of motion for do give the appropriate magnetohydrodynamics. 

2. The canonical quantization of ([18]) leads to the expected current algebra. In particular, 
one finds the equal-time rules 



= irV(x,t)5=^(x-y). (20) 



The charge density, considered as a matrix in the fundamental representation, transforms 
as p — )• ph, h ^ SU{2), p = p^t"". We can thus pick a specific SU{2) transformation g 
which diagonalizes p, 

P = gpdiagg~^ (21) 

so that p'^ = riTr (^ga^g^^t^) . This identifies the dynamical variable ^(x, t) as part of the 
charge density. The eigenvalues of p are gauge- invariant and represented by n. Their flow 
is given by u^. 

For a general gauge group G, the action is given by 

S = -t j'i Tr {qsg-^D^g) - j F(ni, na, . . .) + Sym{A) (22) 

where qg are the diagonal generators of G and jsjsfi = s = 1, 2, . . . , rank{G). 



3 Fluids and gravity 



We now go back to the case of the fluid with no nonabelian internal degrees of freedom. 
As noted before, this case is described by the action ([8]) (or its relativistic version Q). 
Nevertheless, there is something not completely satisfactory about this. The group element 
g belongs to SU{1, 1) and this group has no particular meaning in the relativistic theory. We 
would like to analyze the effect of gravitational or mixed anomalies on the fluid equations. 
The anomalies, as is well known, can be formulated in terms of the Lorentz group which 
acts on the tangent space or in terms of diffeomorphisms. The former point of view requires 
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identifying a Lorentz group action, while the latter can be related to Poincare group action. 
The SU{1, 1) description does not immediately lead to an easily identifiable action of the 
Lorentz or Poincare groups. For this reason, we seek a generalization of the action ([9]); the 
particles underlying the fluid description may or may not have spin. 

3.1 Fluids with spin 

We start by considering an action similar to ^ but with the group element g S SU{1, 1) 
replaced by an element of the Lorentz group, say. A, in some finite dimensional matrix 
representation ^7]. The appropriate action is. 



S[e,u,j,A] = J dete [-if TT{Si2A-^D^{u{e)) A)- F{n)] 

leabcd [ e-Ae^AR'^iio). (23) 



32 7rG 

We have added the Einstein-Hilbert action for gravity as well, written in terms of the frame 

.ab 



field one- form = dx^ and the spin connection lo°'^ = lo^^ dx^ . R""^ is the curvature 
two-form given by 



R"^ = duj"^ + u""" A u:'* = ]^Rf^ dx^ A dx^ 

= \e''^{e-^f^{R^,rpdx^' Adx". (24) 

a;(e) is the torsion free spin connection derived entirely from the metric or equivalently the 
frame fields. It is taken as understood that the contraction of the tangent space indices 
is done with the flat Minkowski metric rj""^. Coordinate indices are contracted, as needed, 
using the metric g^^, = e'^^e^^r]ab, so that = j^fg^u- Further, in (j23|) . is a matrix 
corresponding to the third component of the spatial spin, i.e., equal to the corresponding 
Lorentz generator in the representation corresponding to A. By considering the right trans- 
lations of A by an element of the form exp(iS'i2 ^12)5 we can see that is a covariantly 
conserved current. Under the local Lorentz transformation A, the transformation rules for 
the various quantities are as follows. 

e-^e"^ = A\e\ oo\ ^J\ = A^-^, (A'^)^ - {dAA-^)\, 

R\ = A\R-^{A-y,. (25) 

The variation of the action (j23p with respect to the spin connection uj gives the torsion free 
condition, 

L» A e = 0. (26) 

This can be solved to determine u; as a function of e; we denote the solution as w(e). It 
corresponds to the spin connection derived from the metric via the Christoffel symbols and 
is explicitly given by oj^ = —i uj'^ Sab with 
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In the action (j23p . u occurs only in the last term; in the covariant derivative for A we use 
Co directly, so that 

{A-'D^A)\ = {A-'d^A + A''oj^,Ay b- (28) 

If we had used lv in this term, the condition for vanishing torsion, namely, equation (|26p . 
would be altered. The use of the solution a) is similar to what is done for coupling gravity 
to spin-^ particles, preserving the Riemannian or torsion-free condition. 

In addition to the equation for uj, there are equations of motion for A, and e^. The 
last one corresponds to the field equations for gravity. For the variation of A, we can use 

6{A-^D^A) = -A-^5AA-^Df,A + A-^D^{6A) 

= A-^d^{6AA-^)A + A-i {uf^SAA-^ - ^AA-^w^) A 
= A-\D^e)A. 

where @ = (JAA^^). This leads to the equation of motion, 

^^,.(^/5W) = (29) 

where Q"'' = Tr (5i2 A'^ S''^ A) is the spin density. Notice that the derivative involved in 
this divergence is Levi-Civita covariant and also covariant with respect to the Lorentz group 
action on the tangent space. Similarly, right translations of A in the S'12-direction gives 

-^d^i^jn^V^j^ = 0. (30) 

The equation of motion for j'' becomes 

j^ = -^iTv{Si2A-^D^A). (31) 

The variation of the action with respect to the metric g^^u (or equivalently, the frame 
fields e^) gives the standard terms except for the the variation due to O. The result is 



+ 6S, 



extra 



SS = IJVa Sg^" - ^{R,u - \g,u R) 
T^Q = nF' u^u^ - g^^{nF' - F) 

SSextra = " / ^9 j'' Q^' SO ^,ab (32) 



where we have used (j3ip . The last term can be simplified using 

i^^l^)ab = (e~"^)b (Va'5e/,,a - V^Jea,a) - (e~-^)a (Va(5e^,fe - V^(5ea,b) 

-{e-'Taie-X {Va6e^-V[s6e^re;vmn. (33) 

Here V denotes the derivative covariant with respect to the tangent space and the Levi- 
Civita connection. With partial integrations and using (f29|) . this can be simplified as 



SSextra = j ^ Sg''" V„(j^Q^ + >Q"^) (34) 
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where Q"'^ = Q {e )a(e )^. Thus the energy- momentum tensor is not quite of the 

( f) 

perfect fluid form T^i , rather it is 

T^.u = Tj/J + 2V,(j^Q- + j,g°^). (35) 
The conservation law becomes 

T^f^^^- - 2 {R^p)\ j^Q^^ = (36) 

where we have used pop and identities on the Riemann tensor. The fluid has a spin density 
and the extra term in (|36p is very reminiscent of the coupling of spin and curvature which 
appears at the point-particle level in the Mathisson-Papapetrou equations [TH]. We may 
regard (j36p as the fluid generalization of the latter. 



3.2 Spinless fluids 

As mentioned before, the basic reason for the description given above in terms of the Lorentz 
group is to incorporate easily the gravitational anomalies in the fluid language. Since such 
anomalies, when they occur, are due to fields with spin, we may regard the action (j23p as 
adequate for our needs. Nevertheless, it is interesting at this point to consider an action for a 
fluid of spinless particles so that the energy-momentum tensor has no extra term depending 
on the spin density. Notice also that, from (j3ip . it is the transport of spin which is described 
by the current and not particle number or mass; in other words, we have a fluid of spin 
carriers, spin being their only attribute. It would be useful to consider the flow arising from 
transport of mass. The natural object for this would be the Poincare group, in particular, 
the translations, since spinless particles do have transport of energy and momentum. Since 
we will need matrix representations and traces, it is easier to consider the Poincare group 
as obtained from the de Sitter group 50(4, 1) via a group contraction. In addition to the 
Lorentz generators Sab, we then have Pa = Sa^/R which are the translations (on de Sitter 
space) with 

[p„n] = i^. (37) 

The limit — )• oo corresponds to the group contraction and reduces the algebra to the 
Poincare algebra. A specific matrix realization of the -SO (4, l)-algebra is provided by the 
Dirac 7-matrices and 7a75- 

If g denotes an element of SO {A, 1), then the frame fields for the coset space are given 
by = — iTr(S"*^g~^ dg), and the metric is given by 

ds^ = -TiiS^'^'g-^dg) TiiSa^g-^dg). (38) 

The action for a point-particle is thus 

I[g] = -mj dt^-Tv{S-^g-^g)Tv{Sa5g-'g) 



'-Idt 



2 {Tv{S-^g-^g)f 



m rj 



(39) 
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where, in the second hne, we have used a world-hne metric as an auxihary field. We can 
further reduce this as 



I[g] = Jdt 



J {-iTr{Sa5g 5)) H 7, 



2) 

(40) 



The similarity with the fluid actions we have discussed is evident. This suggest that, for 
spinless fluids, we can use the action 

S = J d^x^ [-ir^Tr(5,5<?-'9^<?) - F{n)] . (41) 

Coupling to gravity is introduced by d^g — D^g = (5^ + 0)^)5, where Cj is the torsion-free 
spin connection as before. The full action is thus 

S = Jdete [-z/^ Tr (^,5 g'' D^g) - F(n)] - -^^^ahcd j e'' A A (42) 

The derivation of the equations of motion will proceed as before. The terms involving g 
will give the energy-momentum tensor of the perfect fluid form, except for the term coming 
from the variation of w; in other words, 

T^, = Tli) + 2 V,(j-g^, + j:Q%) (43) 

where Qc^ = Ti{g Sc5g~^S'^^){e~^)a{e~^)^. Since Sc5 and Sab are orthogonal with the 
trace, Qc^ vanishes unless g Sc5 g~^ generates a term proportional to the Lorentz generator 
Sab- This can only be done via the commutator [5'c5,5'rf5] for terms in g which are of the 
form exp(iS'd5 6"^^ + •••). As a result, Qc^ is of order and vanishes in the contraction 

limit. The energy-momentum tensor then has the perfect fluid form. Thus, the action (j42p 
can describe spinless fluids in terms of the Poincare group realized as the contraction limit 
of the de Sitter group. 



4 Standard Model 

We are now in a position to apply this to the standard model and a fluid or plasma phase of 
the same. For specificity consider the quark-gluon plasma phase for three flavors of quarks, 
u,d,s. In other words, we consider a phase with thermalized u,d,s quarks, so that they 
must be described by fluid variables while the heavier quarks are described by the field 
corresponding to each species. The flavor symmetries, for the moment, will be taken to be 
gauged. We will also neglect the quark masses so that we have the full flavor symmetry 
U{3)l X U{3)r. Thus the group G to be used in 1^ is 

G = SU{3)c X U{3)l X U{3)r (44) 

with individual flows corresponding to the charges. In this discussion our focus is on the 
flavor transport, so we will drop SU{3)c from the equations to follow. Of course, the flavor 
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symmetry is not fully preserved even in the absence of masses due to anomalies. On this 
question, we then have a rerun of the old 't Hooft argument |19j . 

Assume all flavor symmetries are gauged with anomalies canceled by a set of spectator 
fermions. In the fluid phase where u, d, s are replaced by fluid variables, we must then have 
a term in the fluid action which can reproduce the anomalies so that the cancellation with 
spectator fermions still remains valid. (In the usual case where the phase being considered 
is one of confinement and chiral symmetry breaking, this term is the Wess-Zumino term 
constructed in terms of the pseudoscalar meson fields.). In the present case, since we already 
have a description of the fluid in terms of group elements, it is easy enough to adapt the 
usual Wess-Zumino term. Thus our fluid action is given by 



ij^^Tr{gY}D,gL) 



-i k^Tr i^Y^R Df^gnj - i k^Tv ( —g^ D^gn 
-F{n3, ns, msms) + SyMiA) 
+ Twz{AL,AR,gLg\i) 



(45) 



where j^gg apply to U{3)l and A;(^3_g apply to U{3)r and gi G U{3)l, gn G U{3)r. The 
last term is the usual gauged WZ term Twzi^L, Ar, U) given in terms of Ai^, Ar and the 
meson fields U G U{3), and gauged in a way that preserves the vector symmetries, but, for 
our purpose, U is replaced by gidR- Explicitly Twz is given by Witten in [20] as 



WZ 



iN 
' 2407r2 

iN 
487r2 

iN 
487r2 

iN 

"m^ 

iN 
'487r2 

iN 
'487r2 

iN 
'487r2 

iN 
487r2 

iN 
487r2 

iN 
'487r2 



Ti{dUU-^f 
Tr {ALdAL + dAiAi + Al) dUU-^ 
Tr {ARdAR + dARAR + A%) U-'dU 
Tr \{ALdUU-^y - {ARU~^dUy^ 
Tr \Al [dUU-^f + Ar {U'^duf 
Tr {dAidUARU-^ - dARdU^^AiU) 
Tr {ARU-^ALU{U-^dUf - AlU ARU-^{dUU-^f) 
Tr {{dARAR + ARdAR)U-^ALU - {dALAL + ALdAL)U ArU'^) 
j Tr {ALUARU-^ALdUU'^+ARU-^ALUARU~^dU) 

Tr (a]^U~^AlU - AIUArU-' + hj ArIJ-^ AlU ArU'^ Ai \ .(46) 



One of the main results of this paper is that the action given above incorporates all the 
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flavor anomalies in fluid dynamics. Once we have obtained (j46p . we can restrict the gauge 
fields Al^Aji to what is needed for the standard model, namely the SU(2) x U{1) group 
of electroweak interactions. It is straightforward to verify that (06]) does indeed lead to the 
usual chiral magnetic effect. 



5 Currents from anomalies 



5.1 The chiral magnetic and chiral vorticity effects 



For the chiral magnetic effect, we have only a background electromagnetic field turned on, so 
that Al = Ar = —iQ A, where Q is the quark charge matrix given hy Q = diag{^, — |, — i). 
The contribution of the anomaly to the electromagnetic current following from ()46|) has been 
given in |21]. The electromagnetic current for the action (j45p . ()46p is then given by 



j/^ = + ^e''^'"^TT[Q{dMU-^daUU-^dpUU-^) + {U-^d^UU-^daUU-^d,3U)] 



47r^ 



Q'^idpUU-^ + U~^dpU) + -{QdfsUQU-^ - QUQdpU-^] 



(47) 



where is the contribution from the nonanomalous part of the action and we have set 
N = 3. For the main points we want to illustrate, it is sufficient to consider a reduction to 
the SU{2) subgroup; in other words, we will consider basically the up and down quarks. In 
this case, we can take 



U 



V 
1 



(48) 



where y is a 2 x 2 matrix which is an element of SU{2). We take it to be of the form 
V = gL g\i, where gi and gn are now elements of SU{2). The nonanomalous part of the 
current can then be written as 



J3 = - \ ['^3?^3LTr(cr3 5^V3C/L) + ms u|^^Tr(a3C/^ Vsc/ij) 



(49) 



Simplifying ()47|) with the choice of U we have made, 



'J a 



487r 



167r2 

2 + ^Tr (SsL ^3R - 1) 



(50) 



where Xg = gl^dpgi - g'^dpgn and S3L = gl^asgL, = Q^^cF^gR, In simphfying dH]) 
to this form, we have used the fact that there is no rank-3 symmetric invariant tensor for 
SU{2). 

When gi = gn, the last term in ([^U]) . namely, Jg, reduces to 

4 = -£^e^^''''^{dMd^0. (51) 
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This is the chiral magnetic effect discussed in [3]. The quantity VO is related to the fluid 
current for the transport of the U{1)a axial charge. And, correspondingly, in a medium 
in equilibrium, with chiral asymmetry for such charges, we may replace 9 by — I^r), 
where the chemical potentials are for the left and right axial charges. Notice, however, that 
the expression for Jit has added terms when are qr are independent matrices. 



The other terms in (jSOh can be simplified further. First of all, using the Maurer-Cartan 
equations d{g~^dg) + {g~^dg)'^ = 0, we can simplify 



e^^^'^^Trig-^d^gg-^dagg-^dpg) = iTiia^ g-^d,g) [iTT{a^g-^dpg)] . 



(52) 



We can use this to simplify the term e'^''°^Tr(X,^XQX^) in (|50p . Further, from the equation 
of motion for jg and /cg , we find 



iTiiasgl'd^'gL) 
iTT{a3g]^^d''gR) 



-2 "■'^ 



dn3 
dF 



■ u' 



3L 



2 dF 
713 dns 
2 OF 



33 



(53) 



where u'^^ and the flow velocities for the left and right isospin. Using these results, 

the current finally takes the form 



IGvr- 



re^'^"'^a,yl„Tr[(S3L + S3ij)I;3] 



167r2 



s^^'^^^d^Trigl^d^gL QR^d^gR) 



+ 



dF 

5713 



^3L u 



OF 



daUSRB 



(54) 



The last term of this expression involves the vorticity of the flow velocities. This equation 
is thus an expression of the chiral vorticity effect. 



5.2 Mixed gauge-gravity anomaly 

In addition to the flavor anomalies, it is also possible to consider the mixed gauge-gravity 
anomaly in the standard model. The six-form index density which leads to this via the 
descent equations is 

h = i^ {T^F)Tr{RAR) (55) 

where the fleld strength is the one corresponding to the weak hypercharge C/(l)y. The trace 
of the hypercharge vanishes for each generation of quarks by itself, so that this anomaly 
is zero. The possibility of a contribution arises when we consider a plasma where some of 
the quarks, say, the up, down and strange quarks, are in the fluid phase while others, say, 
charm, is to be described by the standard fermion Lagrangian. In this case, for the fluid 
part we would need an effective description. 

There are two choices on how this anomaly can be displayed; we can choose to regard this 
as an anomaly in the hypercharge current or as an anomaly in local Lorentz transformations. 
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For the first point of view, the index density leads, via the descent equations to the effective 
action 



The hypercharge current has the conservation law 



N 



fj,ua/3 



(56) 



(57) 



768 7r2 ^ 

Further, if we choose to regard to as an independent quantity, then the torsion- free condition 
is modified by the Lorentz Chern-Simons term, when this term is added to the Einstein- 
Hilbert action. The more canonical thing to do would be to consider cj in (I56p to be the 
solution w. In this case, with w — )• a), we find for the correction to the energy-momentum 
tensor, 

icorr I92vr2 ^ 

The remaining trace is over the hypercharge values. If we replace 9 by the chemical poten- 
tials, as can be done for the chiral magnetic effect, 

1 

2 



:V> 



Tr(5^0) {R^p 



+ {u ^ a) 



(58) 



Tr(^) 



^(/i^+/zi + /ii) + ^(/x^ + /x|j-2M^; 



(59) 



More generally, we can replace 5^0 by its expression from the equation of motion giving a 
term involving the derivative of the enthalpy function, similar to what was done in (j54p . 
(Since the derivative of the enthalpy at fixed entropy and pressure is the chemical potential, 
this includes the previous case as well.) Thus, depending on the properties of the enthalpy 
function of the fluid, the corrections displayed in ()58p can be nonzero even when = 0. 

The other possibility is to consider the index density as leading to anomalies in local 
Lorentz transformations. We can use an element of the Lorentz group, identified as the 
fluid variable A of section 3, to write the Wess-Zumino term. The transformation of fields 
of the relevant flelds is given by 



A^A3 = gA 
^-dgg-\ R 



e ^ e" = ge, 
UJ ^ Lo^ = guj g 

The Wess-Zumino term may then be written as 

N 



R3 



9^9 



(60) 



^wz 



1927r2 

N 



Tr(F) 



T:i[ojdoj + -J^ 

O 



Tr ( ndVL + -^^ 



I Tr(F)Tr(dAA-i^) + iTr(F)Tr(dAA-i)3 



(61) 



where Cl = A^^ [dA + coA]. Once again, if we regard u as independent, then this leads to 
a nonzero torsion proportional to the spin-density. The equation of motion for ui, starting 
from (p3|) and adding (j6T]) . can be reduced to the form 

{T^,r = {M»r'ie-')l-{M('r'{e-')l + {M'r''{e-'me-'fde-^ 
N 



cd 



192 TT^ 



^f^ua^^abcd Tr(F^^) Tr(5„A A'^ S, 



ab 



(62) 
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where (Tf^u)"" = {D^CyY ~ {D^e^Y is the torsion tensor. 



In the case when we use uj in place of a; in (j6ip , we get corrections to the equation of 
motion. For variations corresponding to the right translations of A by a term proportional 
to 512, we find 

— D^i^f) Tt{SI^) = --^ Tt[F,,) 5,Tr(5i2 A'^ DpA). (63) 

The equation for the left translations of A by an arbitrary infinitesimal Lorentz transfor- 
mation is 

— [^f A5i2A-i] = TV(F^,) [R^p - D^{DpKK-^)] . (64) 

(Of course, the two equations, ([63]) and (fM|l . are not completely independent.) 

The variation of ()6ip with respect to the frame field will yield the correction to the 
energy-momentum tensor. This is given by 

nleor-. = -^^eA,aV;3[Tr(F^,)(Tr(5,AA-iS^'^)e'^'^-^ + Tr(a„AA-iS'5^)e'^'^"'^ 

1 iN 
96^ 



— 7ir2^A,aV;3 Tr(F^,)TV(a„AA-i5^-)e^^"^ (65) 



where S°^^ = 5°^(e~^)"(e~^)'^. The terms in the first line of this equation leads to a 
symmetric energy-momentum tensor when written in terms of the coordinate components, 
by multiplying with {e~^Y°' . The term in the second line leads to an antisymmetric term. 
This is to be expected. We know that a symmetric energy-momentum tensor is necessary 
for the conservation of the current corresponding to the Lorentz transformations. In the 
present case, the anomaly implies that this current is not conserved. The antisymmetric 
term is a manifestation of this property. Since the Einstein tensor R^u — ^g^ivR is symmetric, 
this leads to a problem with the Einstein equations. The proper way to understand this is 
to realize that there is another term in T^y due to the quarks we have neglected, say the 
charm quark in the example we have been using. Since the latter field by itself also leads to 
a mixed anomaly, the full energy momentum tensor which is the sum of T^L/ , T^j^ from (j65p 
and T^i, from charm will together be symmetric, the anomaly part from the charm quark 
canceling the antisymmetric piece of (j65p . 



6 Discussion 



We have obtained a very general formalism for incorporating the effects of anomalies in 
hydrodynamics. As mentioned in the introduction, being a formalism based on symmetries 
rather than calculations specific to any particular assembly of material particles, this is quite 
general and is expected to be valid beyond weak coupling or near equilibrium conditions. 
The specific choice of the fluid will be reflected in the choice of the enthalpy functions. 
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A few clarifying comments are in order. It is important to realise that in any fluid where 
the particles which constitute it carry a variety of quantum numbers, we can have a number 
of different flow velocities. This is evident from the action (j22p where we have flow velocities 
for each diagonal generator of the group. This point seems not to be adequately emphasized 
in the literature. It is also useful to visualize this as follows. Consider two quarks and two 
antiquarks in a fluid. We could have them forming a color singlet and moving in the same 
direction. This gives a mass/energy flow but no color flow. We could visualize a, qq pair 
forming an octet state and moving together in a certain direction while the other q q pair 
form a singlet. This gives a nonzero color transport rate different from the mass/energy 
flow. We could also envisage subsets of particles forming different spin states giving a spin 
flow velocity, possibly different from the mass and color flows. When we consider massless 
quarks, the L, R quantum numbers are independent quantum numbers with independent 
velocities possible. 

Specifically for the flavor part, we can have independent u^l and u^r. These need not 
coincide even when gi = gn for two reasons: The local charge representation is determined 
by n3, ms and these need not be the same even when gi = Or- Secondly, the enthalpies can 
be different as well. 

If a calculation is carried out in a specific medium, the results obtained would be for the 
appropriate enthalpy function. For example, if we take a massless field, then the relation 
between pressure and energy density corresponds to an enthalpy function F ~ n^^^ . In this 
case, (-F')^ is of the form tn?/'^ . Notice that for the vorticity term in (j54p there is a prefactor 
proportional to this. If, in addition, we take n ~ T^, as is appropriate for a relativistic 
gas, then the prefactor gives a term. This may give a point of correspondence with 
the results in p^. However, we should expect a contribution even at zero temperature, 
since the structure of the anomaly has to be reproduced correctly in the fluid language; this 
is evident from section 5.2. The derivative of the enthalpy function is also related to the 
chemical potentials, when the latter is introduced. Our formula (j54p is thus similar to the 
results in [22J as well. 

Regarding the use of the Wess-Zumino term for anomalies, the specific choice of T\yz 
specifies the nature of the currents being discussed. (This point is moot for our discussion 
in sections 2 and 3, since we have not introduced anomalies yet.) We have used the form 
(|46p which gives expressions invariant under the nonanomalous vector gauge symmetries. 

The Wess-Zumino term was also used to obtain anomalies for chiral superfluids in [23], 
although the formalism is very different from ours and the emphasis was on baryonic and 
axial currents. (This article came to our attention after this paper was completed. We 
thank the author for correspondence on this.) However, we may note that equation (58) of 
[23] is similar to our ([5^ if our u^l and u^r are related to the different superfluid velocities 
introduced in that paper. 

This work was supported by U.S. National Science Foundation grant PHY-0855515 and 
by a PSC-CUNY grant. 
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